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Abstract 



> 

' We investigate the interaction of an atom with a multi-channel squeezed 

. vacuum. It turns out that the light coming out in a particular channel 

I can have anomalous spectral properties, among them asymmetry of the 

' spectrum, absence of the central peak as well as central hole burning for 

particular parameters. As an example plane-wave squeezing is considered. 
In this case the above phenomena can occur for the light spectra in certain 
directions. In the total spectrum these phenomena are washed out. 



C ■ PACS numbers: 42.50.Dv, 32.80.-t 

1 Introduction 

^ I The interaction of squeezed light with atoms has found considerable theoretical 

I attention in the last years. As pointed out in the seminal paper of Gardiner 

[Q, squeezed light, either in the white noise limit or broadband colored noise 
[H, can drastically alter the radiative properties of atoms and can, in principle, 
reduce substantially the spectral linewidth of emitted light. Gardiner's paper 
has stimulated a large amount of work on the interaction of squeezed light with 
atoms. We refer to the detailed review of Parkins P] for references and an account 
of work up to 1993. For more recent work see, e.g., Refs. Q and references 
therein. 

Gardiner ^ considered a one- dimensional like situation in which only a 
single channel is squeezed and where the atom couples only to this channel. Here 
a channel means a set of modes ranging over all frequencies with the remaining 
quantum numbers characterizing the modes kept fixed For example, a channel 
may consist of all modes with fixed angular momentum and parity for a multipole 
expansion of the field, or of all plane wave modes with a given direction and 
polarization. 
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An investigation of resonance fluorescence spectra of atoms interacting with 
a (single-channel) squeezed vacuum plus a laser was initiated by Carmichael, 
Lane, and Walls (for further references see the review and Refs. 0, 0). 
Smart and Swain M, 3 have pointed out the existence of interesting structures in 



these spectra. In Ref. [|Tl|, |T0[ multi-channel squeezing and associated correlation 
functions for three-level atoms were studied. 

In this paper we consider the spectral effects of a multi-channel squeezed 
vacuum in the white noise limit on a two-level atom. For the atomic correlation 
functions and the total spectrum of all outcoming light a multi-channel squeezed 
vacuum leads to analogous results as a single-channel squeezed vacuum with 
appropriate parameters. However, in a multi-channel situation one can observe 
not only the total spectrum but also the light spectrum in individual channels. It 
turns out that — due to interference of the (quantized) light scattered from the 
atom with the squeezed vacuum — these spectra can show unexpected features 
which are not visible in the total spectrum, e.g., a possible asymmetry, absence 
of the central peak as well as central hole burning for particular parameters. By 
the same arguments as in Ref. we expect these features to persist also for only 
approximate white-noise squeezing. 

In Section ^ below the general case of multi-channel squeezed white noise 
interacting with a two-level atom is treated. The spectrum is calculated in terms 
of a background, scattered, and interference part. 

In Section |^ we treat in detail an example in which the channels consist 
of plane waves with fixed directions and polarizations. In this case one has a 
divergence of the background term when calculated in terms of photon numbers, 
and we therefore use the spectral Poynting vector to calculate the spectrum for 
given position and direction of observation of the spectral analyzer. In this case 
the above phenomena like asymmetry, central hole burning etc. can occur for the 
light spectra in certain directions. 

In Section | we discuss our results in detail, in particular the question of 
interference, and point out a possible connection, probably more formal than 
directly physical, with the results of Ref. H]. 



2 Spectra for general mult i- channel squeezing 

We consider a two-level atom coupled to the electromagnetic field in three- 
dimensional space. The free Hamiltonians of the atom and the field are given 
by 

poo 

Hx = hujQa^a^ Hp = / dujy^ hu al^{uj)aa{uj) , (1) 

where ujq denotes the atomic transition frequency, 0"+ = |+)( — | and cr^ = |— )(+| 
are the atomic raising and lowering operators and the aa{uj) are the annihilation 
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operators of the field obeying the commutation relations 

[aai{uJi),ai^{uJ2)] = Sa^ai^i^^l -^2) ■ (2) 

The index a stands for all quantum numbers of the chosen modes apart from 
their frequency and thus characterizes a channel. These are, for example, parity 
and angular momentum quantum numbers if the multipole expansion of the field 
is used, or the direction of propagation together with the polarization for plane 
wave modes. In the rotating-wave approximation, the interaction Hamiltonian 
may be written as 

roo 

HAF = —f= / duj^g^{uj)a)^{uj)o' + h.c. (3) 
V27r JO a 

with possibly complex coupling coefficients gaij^), which we decompose as 



g^{uj) ^ ^^^{u) e^^'^(-) . (4) 
Going over to the interaction picture leads to the Hamiltonian 

H^{t) =inY,^ [bUt)(r- - b^{t)a+] (5) 

a 

with 7q, = 7a(co'o) and 



1/2 



271 Jo XlaiuJo) J 

Assuming that the requirements for applying the Markov approximation are sat- 
isfied we replace the second factor in the commutator 

Ms), blit)] = S^p-4 r du e-^(---°)(^-*) ^ ■ 6{s - t) (7) 

Zir Jo la{i^0) 

by a 5-function (cf. 

The radiation field is supposed to be initially in a pure broadband squeezed 
vacuum state with the atomic frequency as central frequency. In general, for 
such a state the second order moments of the ba{oo) read in the white noise limit 
0,i 



( &a(s)&/3(t) ) = M^fs 6{s - t) { bi{s)bf,{t) ) = No^p 5{s - t) (8) 

with 

M^ = M U^=U M^M=N{U+l). (9) 
If one omits factors of the form [70,(0;) (t^o)]^''^ this leads to 

(a,(cui)a^M) = A<a/3 5(2cuo-a;i-cu2)e'['^"(-^)+'^''(-^)l 
{aiMapiu^)) = Aro^5(a;i-cu2)ei[-<^'^(-^)+<^^(-^)l (10) 
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for the moments of the aa{uj). In the following Eq. (j^) will be used as definition 
of squeezed white noise. 

We shall assume in the following that Ai and A/" are diagonal for the given 
modes, 

Maf3 = Sap M„ AfafS = 5^(3 • (11) 

One can imagine a state of this kind as being produced by independently squeez- 
ing modes with different a, e.g., by coupling them to different parametric oscilla- 



tors [jT3[. For a pure state, as considered here, the assumption ([TI|) seems not to 



be very restrictive. See [|T4|, where for a finite number of modes the question 
is discussed to what extent second order moments can be simplified by choosing 
appropriate modes. 

The fluorescence spectrum of an atom illuminated by squeezed white noise 
has been studied by Gardiner |]l| and has since become textbook material [^]. In 
Gardiner's article and in a large part of the following work it was supposed that 
initially only a single channel of the radiation field is squeezed and that the atom 
couples only to this channel. In our notation, this means there is a particular a, 
a = 0, say, with 

Ma = M5a,o Na = NSa,o la = lSa,o. (12) 

Such a situation will, a little imprecisely, be called one- dimensional in the follow- 
ing. 

We now want to calculate the spectrum Sa{uj) that would be observed by 
a spectral analyzer coupled to the modes {(y,uj) with a fixed. (For plane wave 
modes this would simply be the spectrum of light with a certain polarization 
observed in a certain direction.) One could of course also observe and determine 
other spectra, e.g., the spectrum in a channel which is a superposition of different 
a's. Our procedure easily carries over to this situation. To determine Sa{uj) we 
shall adapt the procedure of Gardiner |^ to multi-channel squeezing. 

At a finite time the spectrum is proportional to the expectation value of the 
photon number operator al^{uj)aa{uj). The stationary spectrum is given by 

Saiu)=\im -i- Tdti rdt2e-i('^--°)(*^-*^)u;„(ti,t2) , (13) 
T^oo ZttI Jo jo 

with the two-time correlation function 

Wa{tl,t2) = { &l(tl)out ba{t2)ont ) , (14) 

and can be written as the Fourier transform 

S^{uj) = — r dre-^("-"°)"w„(r) = - Re /"^ dr e-'("-'^«)"w^(r) (15) 
2tx J-oo tt Jo 

of the stationary correlation function 

Wa{r) = lim Wa{r + t2, ^2) ■ (16) 

12^00 
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The operators &a(s)out appearing here denote the hmits 

ba{s)out = lim ba{s)t , (17) 

t — >oo 

where the subscript t stands for the time evolution in the interaction picture, 

Xt = Uj{t,OyXUi{t,0) . (18) 
Since the ba{s)f obey the equations of motion 

^^Ush = \ [Hi{t)tMs)t\ = 5it-s)^a:^ (19) 

one has 

bais)t = ba{s)+i3it~s)^as , (20) 
where denotes the Heaviside function, and in particular 

bais)out = bai-S) + ag . (21) 

Inserting this into Eq. (^3|) yields a decomposition of the spectrum into three 
parts, corresponding to the correlation functions 

w\{hM) = v/7^(&l(ti)a^- +4^6„(t2)) , (22) 

which will be called the background, scattered, and interference part, respectively. 
For the background part one finds immediately 

wlihM) = wlih - t2) = NJ^Jiti - t2) 27r5^(^) = NJ^^ , (23) 

as expected for the spectrum of white noise. Although 6aa = 1, it has explicitly 
been kept since it becomes divergent if the index a is not purely discrete, e.g., 
for plane wave modes. For such modes, a state with Ai, Af being diagonal is 
an idealization just as a plane wave coherent state. It will be shown in the next 
section how this problem can be bypassed in a simple physical way. 

The 'mixed' correlation function w^(ti,t2) can be reduced to an expression 
containing only the atomic operators by means of the formulae 

{bi{ti)%) = m-h).^{N4al,a:^J-M:[a:^^,a^J) (24) 
(a+6,(t2)) = ^9(ti-t2)-v^(iV4a+,a^-]-M„[a+,a+]) (25) 

which can be derived in a similar way as in the one- dimensional case. 

The remaining correlation functions of the cr^ can be evaluated in the station- 
ary limit using the quantum regression theorem and the atomic master equation. 
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Since the calculations for the one-dimensional case in can be carried over, 
except for some obvious modifications, we just give the results. The form of the 
master equation 

p = ^(iV + l)(a-pa+-^{a+a-,p})+7iV(a+pa"-i{a-a+,p}) 

-7M a+pa+ - 7M* a'pa~ (26) 

remains unchanged, but the parameters 7, M, and now are defined as 

7 = E7a M = 5:^M„ N = Y.^-^N^. (27) 

The possible values of M and are restricted by the inequality 

\M\^ <N{N + 1) , (28) 

which follows from the relations |Mc,p = No,{N^ + 1) [cf. Eq. (D]. 

In the one-dimensional case 7, M, and N defined above agree with those 
appearing in Eq. (|T^) , and for pure states — which we have been considering 
above — the equality sign would hold in Eq. (p8|) . The full range of the parameters 
M and N can also be realized in this case if one uses mixed states. 

By absorbing a phase into the atomic states if necessary, M can be chosen 
real and positive. With this convention one gets for r > (the relation w^'^(r) = 
w^'^(— r)* yields the corresponding expressions for r < 0) 



la N 

2 2N + 



I / \ 7q 

wJt) = 



2 2N + 



-{e-^+^ + e-^-^} (29) 
^ {e-^+^(iV„ + M„) + e-^-^(iV„ - M„)} , (30) 



where 



7± = (iV±M + l)7. (31) 



The complete spectrum is the Fourier transform of 



2 2A^ + 1 

X [e-^+^{N -N^- M„) + e-^-^{N - A^^ + M,)} , (32) 



7± 



2N + 1 ^ {uj- cJo)2 + 7± 



+ 2M7„72 



U — UJq 



[{uj - wo)^ + 7+] [(^^ - '^o)^ + 7 



(AT-AT^^ReM^) 

ImM„. (33) 
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For the spectrum of all modes, S{uj) = J2a<Sa{^), one obtains 

2.S(.) ^ A' + ^ + - (.-J„K,J ■ (3^) 

an expression that coincides with the spectrum in the one- dimensional case. This 
shows that, as long as only this kind of spectrum is observed, all states satisfying 
( p!TD are equivalent to mixed states of the one- dimensional type (|I2|); c.f. the 
remark following Eq. (^). 

In contrast to S{uj), Sa{uj) shows new features, namely in general the spectra 
are asymmetric and the sign and the relative weight of the peaks depend on Ma 
and Na- The asymmetries are caused by the phase of Mq, (this phase has a 
physical meaning as it is the relative phase of Mq, and M) since the symmetric 
and antisymmetric parts of Sa{ujQ + uj') are the Fourier transforms of the real and 
imaginary part of Wa{T), respectively. 

In figure the spectrum (|3^) is plotted for Na = N and increasing values of 
(p = arg Ma- For (p = one obtains the spectrum of the one-dimensional case 
consisting of two Lorentzians, a positive narrow peak of width 7_ and a negative 
broad peak of width 7+. For < (f < n the spectrum is asymmetric with respect 
to ujq, the asymmetry being maximal for (f = 7i/2. The relative weights of the 
symmetric contributions to the spectrum decrease and vanish for (f = tt/2. For 
n/2 < ip < IT, the weights increase again but the signs of the two Lorentzians 
are interchanged. The spectra for vr < (y9 < 27r coincide with those for 2tt — cp 
mirrored dX uj = ujq. 

For Na = (and consequently Ma = 0) the background and interference 
parts of the spectra vanish. As the scattered part does not depend an Na and 
Ma, the spectra take the shape of the dotted line in figure |l] (the dashed line now 
representing S = 0). 

Spectra for Na > N are plotted in figure |[ As compared to those with 
Na = N, the positive peaks are attenuated while the depth of the negative peaks 
increases. 

3 Example: Squeezed plane waves 

To be more specific, we shall now deal with plane wave modes in greater detail. 
These modes are specified by their wave- vector k and their polarization A = 1, 2, 
i.e., hj uj = ck and a = (k. A) with k = k/A;. In the dipole and rotating-wave 
approximation the interaction Hamiltonian is given by 

i^AF = -d • E(-) a- - d* • EW a+ , (35) 
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where E'^+^ = E*- -"^ denotes the positive frequency part of the electric field at the 
position r = of the atom, 

(27r)2 J \Zeo/ ^=1,2 ^ 

and d = (— | D |+) is a matrix element of the atomic dipole operator D (by parity, 
the static dipole moments in the states |±) vanish). 

In order to apply the formulae of the previous section, the sums over a and 
the Kronecker symbols Saia2 have to be replaced according to 

(d^-k 6x,xj\kl,k2) , (37) 

where d^k denotes the area element on the unit sphere and (5^(ki,k2) the corre- 
sponding 5-function. By means of the identity 

S\k, - ka) = A;f 2 S{k, - 5\ki, k^) (38) 

one can write 



^3 

[aki,Ai, aj,^ = (5^(ki - k2)6x^\2 = — ■ S^{ki, k2)6xi\2 ■ S{uji - UJ2) (39) 



and by comparing with Eq. (^ one sees that operators with commutation rela- 
tions analogous to those of the aa{uj) can be defined by 

^k,A(^) = ak,A • (40) 

The coupling coefficients appearing in a representation of i^AF in the form of 
Eq. (H) are therefore given by 



»» = ifd:^) .tA.d.,/^e>*M. (41) 



Note that their phases do not depend on u. 

As already noted, a squeezed state with Ai, M being diagonal for plane 
waves leads to a divergent expression in Eq. (^). One could try to circumvent 
this problem by subtracting the divergent background part, but a physically more 
satisfactory solution is to use an improved definition of the spectrum which avoids 
divergences automatically. 

A suitable quantity for modeling the spectrum actually observed in an exper- 
iment is the spectrally resolved energy fiux through a (small) surface A centered 
at r = ro, 

SA,roi^)=[ do-.(S(r,a;)), (42) 

•J {A, ro) 
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where the operator S(r, u) represents the 'spectral Poynting vector', i.e., the spec- 
tral energy flux density. More realistically, one could also use a direction-sensitive 
spectral analyzer, e.g., an analyzer admitting only directions in a certain cone in 
k space, such that only radiation with directions from this cone is observed. It 
can be shown that, as physically expected, the scattered and interference parts 
of the spectrum are not influenced by this as long as the cone contains the line 
between the atom and the detector. For the background part the directional 
selection would simply result in a restriction of the integration over k in Eq. (^) 
below to this cone. 

For the operator S(r,co') we shall use the expression (cf. |T^) 

S{r,uj) = lim ^E5r)(r,cu) X Bijt\r,uj) + h.c. , (43) 

T — >oo 2 

where Ej^''(r, cj) — and analogously 'B^\r,uj) — is defined by 

'2tx h 



E^.^) (r, c.) = ^ dt e±-* E(±) (r, t) , (44) 



E(='=)(r,t) and B'^='')(r, t) being the Heisenberg operators of the positive and neg- 
ative frequency parts of the (transversal) electric and magnetic fields, 



(27r)i J V2eo/ 



A=l,2 
1/2 



BW(r,t) = -^/d^k(^)' E^kx4,e''^-a,,.(t). (46) 



A=l,2 



Introducing the correlation function [|Tj] 

WA,nitxM) = ^oC^ [ d^T.e-i-o(*^-*^)(E(-)(r,ti) xBW(r,t2) 

J {A, ro) 

-B(-)(r,ti) xE(+)(r,t2)) (47) 

and its stationary limit WA.,ro{'T) = limt2^oo w^,ro(''" + hyh), <SA,ro{^) can be 
written as a Fourier transform in the same way as Sa{uj) in Eqs. (|T3|) and (|T^), 

SA,ro{oo) = hm -i- Tdti rdt2e-(---°)(*-*^)^^,ro(tl,t2) 

T^oo Znl Jo Jo 
1 

= -Re/ dre-'("-'^")"w^roM • (48) 

71 Jo 

Since the time evolution operator U{t) = exp (^—ih^^{Hx + Hp)t^ Ui{t,0) 
transforms ba{s) into 

f/(t)t 6,(s) U{t) = e-^^o* b^{s + t)t , (49) 
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it follows from Eqs. (|^) and PT| ) that the a^^xit) are related to the b^^(t + s)f by 

ei^ii,. e--o(*+^) bi^,{t + s)t = ^ Tdu (-Yuc-I e"'-^ a^,,{t) (50) 

\/2n Jo \uJo/ 

with p = 3/2. Within the scope of the Markov approximation, this equation 
remains valid also for p 7^ 3/2 since the factor uj/uq can be replaced by unity just 
as 7a(^)/7aK) in Eq. ® (cf. Ch. 8.1). 

Because of the linearity of the above relation, the decomposition of the x{s)f 
according to Eq. (|20|) , 



^ h,xit + s)t^''>+bi,^,{t + s)t^'\ (51) 
leads to an equivalent decomposition of the ak,A 

(t) as well as of the fields E(+)(r, t) 
and B^+)(r, t). Thus the correlation function Wy^,ro(^i5 ^2) can be written as a sum 
u!A,ro = "^S.ro + '^Aro + ""^^ro aualogous to Eq. (^) for Wa- In the following, we 
shall denote by ( . . . )^'^'^ the sum of all the terms in (...) that contribute to 

B.I.S • , 

^Xro ' i-e-) we set 

{b{b,)^ = {bY%2^'^) {b\b,r = {b{^'V^), 

{b\b,y={b{^%,('^ + b{^%,('^) , (52) 

with bi standing for ^^(tj + Si)f.. 

Using Eq. ( |50D the correlation function of the background part can easily be 
calculated. With 

2 H(ki, Ai; k2, A2) = X (k2 X e^ J - (ki x ,J X 4^ (53) 
one finds 

WAroiti-t2) = ^ / do-. / d^ki / d^k^ S(ki,Ai;k2,A2) 

X exp (-i0k^,Ai + %2,A2) (i^oC~^(ki - kg) • r) 

= ^(^1-^2)-^/^ d^. /d^k^iVt^.k, (54) 

Aq JiA,ro) J ;^ 

where Aq = 2'n'c/uJo is the wavelength corresponding to the atomic transition 
frequency. 

The calculation of the interference and scattered part can be performed with 
the aid of the asymptotic expansion [l^, 

Id^ke-^-^lk) ™ ^{e-i^^/(f)-e''=^/(-f)} + 0(r-2) (55) 



10 



being valid as long as f(k) is sufficiently smooth. Note that this condition is 
violated for the background part pO| calculated above where /(k) oc 5^(k, k'). 
By ffist applying (^) to either integral over k in w^^ro' keeping only terms of the 
leading order r~^, and then using Eq. (|50|) , one obtains 

X exp (-i05if,Ai + i06f,A2) exp (iuJoc~^{^i - 6)^^) 

X ( blrM (^1 - %f ,A. (^2 - c~%r)t^ ^ (56) 

Due to 

H(k,Ai;-k,A2) =0 H(k,Ai;k,A2) =k5AiA2 (57) 

the sums are actually running only over = ^2 and Ai = A2. Further, since 
(. . .)^'^ contains at least one factor {}{ti — + c~^^jr) = ^{^i), only the term 
with ^1 = ^2 = +1 survives. For a surface which is sufficiently flat and whose 
diameter is small compared to its distance to the atom, r can be approximated 
by ro, the integration over {A, tq) resulting in a multiplication with the oriented 
area A. = ^0) finally have, with w^f defined as in Eqs. (pOj) and 

(ID, ' 

^Aro(^i'^2) = ^^0 ^2^^° E^^o^l^i -c'^ro, t2 - C^Vq) (58) 

and 



^0 A 



By inserting Eq. ( Pl) one explicitly finds for 7, M, and in Eq. (|27| ) 



/ )3 lrl|2 

^ = 1^. (60) 



(this is the Einstein A coefficient for a dipole transition) 

A 

and the three parts of the spectrum read 
2n 



M 

Stt 



/ d^kE Isfc,, . d|X, AT = A /d^kE Isfc,, . d^iV, , (61) 



27r s 3A-fo A^7 /, - |2\ 7± .^^n 

27r I 3 A • fo v-^ I ^ 
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X • 



with Ao the wavelength of the atomic transition frequency. We note the direction 
dependence of the spectrum and the dependence of iS^^^o '^Aro- 

Eq. (^) shows that, as far as the interference and scattered parts are con- 
cerned, the spectrum in Eqs. (^) - ( |5^ is a sum of spectra for fixed polarization. 
The latter correspond to a spectrum for fixed a as in Section ^ and can be ob- 
tained by omitting the sum over A in Eq. (|6^) and replacing the factor 1 — |d • foP 
by \ero,\ • d| in Eq. (||). 

The basic features of the spectra are the same as in the general case of Sec- 
tion B. This will now be discussed in more detail. 



4 Discussion 

The spectra in Eq. (pSf), as well as those in Eqs. (|62|) - (@), are a sum of a 
background term S^, a scattered part and a term S^, which can be identified 
as an interference part. That this term is indeed due to interference (note that 
all radiation is quantized) can be seen in various ways. Formally, this is already 
suggested by Eqs. (^) and (p^. The last term in Eq. ( [21| ) is due to the presence 
of the atom since it vanishes for 7q, = and Eq. (^) is due to the cross terms of 
this with 6„[PT . 



The interference nature of becomes yet more transparent if one uses the 
spectral Poynting vector for the calculation of the spectrum, as in Section ^. As 
pointed out there, the results for and remain the same for a direction- 
sensitive spectral analyzer, as long as the analyzer points in the direction of the 
atom. If the direction-sensitive analyzer does not point in this direction one 
shows by the same arguments as in Section ^ that and become zero. This 
is physically expected since it means that the corresponding light originates at 
the site of the atom. Furthermore, if the vacuum is squeezed only for directions 
in some cone C, i.e., if N^^x vanish for k not in C, then vanishes if 

the direction fo from atom to analyzer is not in the squeezing cone C, since then 
A^f„_A and Mf^ a in Eq- (@) vanish. Moreover, depends on N^^^ and M^^^ only 
for k = f [except for the general dependence of N and M in Eq. (0)], and this 
can be interpreted as the fact that scattered light interferes only with incident 
light traveling in the same direction, just as for classical light scattering. This 
can be traced back to Eq. (|55D . 

The spectra in Eqs. ( p^ ) - (Q) depend on the direction f o of the analyzer. As 
seen from Eq. (|63|), outside the squeezing cone the spectrum consists of a narrow 
Lorentzian sitting on top of a broad Lorentzian. This is a special case of the 
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dependence on a in Eq. (pS]). As remarked earlier there is no interference outside 
the squeezing cone. The analogous fact is true in the general case, as seen in 
Eq. (|30D for S\ which vanishes when the channel a is not squeezed. 

For multi-channel squeezing one has more parameters available than in the 
single-channel case. In addition to and M of the single-channel case one 
now also has A^^^ and (or A^g-A ^'^'^ ^icA plane- wave model), with 

|M„p = Na{Na + 1). The range of M and A^ is restricted by the inequality 



\M\'^ <N{N + 1), (65) 

where the equality sign holds, as easily shown, if and only if A^^ = A^ and = M 
for all a with 7^ 7^ 0. Only in this case does the spectrum Sa{uj) in Eq. (^) 
contain a peak which becomes increasingly narrow for increasing A^ as seen from 
the definition of 7_ in Eq. ( pTD which for constant ratio m = M[N(N + 1)]^^/^ 
can be written as 



7- 



iN + l)il-m) + ^ + 0{N 
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For m = 1 this not only recovers the one- dimensional like case of Gardiner [|l| 
with its interesting narrow peak ||2^ , but also shows that a subnatural linewidth 
and the above new features are mutually exclusive. A similar situation is found 
in Ref. for single-channel squeezing with an additional laser. 

If the equality does not hold in Eq. (p5|), then 7_ increases with increasing 
A^ and the corresponding peak cannot become arbitrarily narrow. But as long as 
MV[A^(A^+ 1)] is not to small (the spectra in figure |I| and belong to MV[A^(A^+ 
1)] = 0.75), there are other interesting features in this case. First of all, if 
is real then the spectrum is symmetric. However, if is complex (since we 
have chosen M as positive this actually means Ma/M complex) then ImM^ 7^ 
and the spectrum is asymmetric. But even if Mq, is real, new phenomena occur, 
as seen in figure ^ and ^ for = 0. In this case the negative contribution 
can substantially exceed the positive one in absolute value, as shown in figure ^ 
{(p = 0), and for negative Ma/M the central maximum can be absent completely. 
The same is true in Eqs. (|63|) and (^4[) for the plane wave model. 

If Mq/M is complex then the last term in Eq. (^), which comes from the 
interference part, makes the spectrum asymmetric. For the plane- wave model 
this is seen in Eq. (0). This asymmetry is a new phenomenon which does not 
occur in the one- dimensional case. Various typical spectra are shown in figure |1] 
and 1^. 

It seems that the relative phases of M and Mq, play a similar role for the 
spectra as the introduction of an additional laser with a relative phase in the one- 
channel case, which has been investigated in Refs. The spectra obtained 
in Ref. |^ resemble those in our figures. The similarity is particularly striking 
for figure [53]. Spectra with a 'pimple' and central hole burning occur also 
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in Ref. 0, and we have a similar sensitivity of this effect on the parameters 
(in figure ^ the 'pimple' disappears for ^ = 0, corresponding to the cancellation 
of the contributions from the scattered and interference part, respectively). It 
seems to us, however, that this is more a formal mathematical similarity of two 
physically distinct situations since in both cases one has sums and differences 
of Lorentzians and the possibility to adjust various parameters. Physically, the 
spectra calculated in Ref. p belong solely to the scattered light from the atom, 
which is driven by the combined field of the laser and the squeezed vacuum, and 
therefore these spectra correspond to detection directions away from the driving 
fields. In our case, however, these spectra result from the interference of the 
radiation emitted by the atom with the squeezed light traveling away from the 
atom, and in these directions the combined light is spectrally analyzed. 

In summary, we have shown that in the case of an atom in a more general 
squeezed vacuum the spectrum can show new phenomena compared to the one- 
dimensional like case, among them absence of the positive peak and asymmetry. 
In the plane-wave model the shape of the spectrum can become direction depen- 
dent which, in the general case, is translated into a dependence. This dependence 
is due to an interference effect of the squeezed light with the scattered light. 
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Figure 1: Typical spectra 5«(tu) for A^-^ = = 0.25, = 0.75 A^(A^ + 1) and 
different pliases ip of Ma (solid line). The dashed line denotes the background 
part of the spectrum, the dotted line the sum of the background and scattered 
part. The spectra belong to the largest values of 70/7 compatible with Eq. (pTf) . 
The y axes are in units of Na/2T[. 




Figure 2: The spectra of figure [1|, but with A^^. = 8A^, = 0.25, and = 
0.75A^(A^ + 1). 
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